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Abstract 

We theoretically study the tunnelling current of a single electron transistor (SET) 
under optical pumping. It found that holes in the quantum dot(QD) created by optical 
pumping lead to new channels for the electrons tunnelling from emitter to collector. 
As a consequence, an electron can tunnel through the QD via additional channels, 
characterized by the exciton, trion and biexciton states. The binding energy of exciton 
complexes can be determined by the Coulomb oscillatory tunnelling current. 
PACS numbers: 73.63.Kv, 73.23.Hk and 78.67.Hc 
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Recently, the spontaneous emission spectrum of a single quantum dot (QD) has been suggested 
as a single-photon source, which is important in the application of quantum cryptography 1-4 . Ex- 
perimentally, the spontaneous emission spectrum typically exhibits coexisting sharp emission peaks, 
which have been attributed to the electron-hole recombination in the exciton, trion, and biexciton 
formed in the QD 5 . Nevertheless, it is difficult to experimentally determine the binding energy of 
exciton complexes due to lack of free electron- hole recombination in the QD. In this letter we propose 
that the tunnelling current of single electron transistor (SET) under optical pumping can be employed 
to determine the binding energy of exciton. Our calculation is based on the Keldysh Green's function 
approach within the Anderson model for a two-level system. We find that the optical excitation cre- 
ates holes in the QD, which provide new channels (via the electron-hole interaction) for the electron 
to tunnel from the emitter to the collector. Consequently, an electron can tunnel through the QD via 
four additional channels, characterized by the exciton, positive trion, negative trion, and biexciton 
states. Each addition channel can generate a new oscillatory peak in the tunnelling current charac- 
teristics in addition to the typical peaks caused by the electron-electron Coulomb interactions. The 
binding energy of exciton complexes as well as electron charging energy can be determined by using 
the tunnelling current as functions of gate voltage. 

The system under the current study considers a single quantum dot (QD) sandwiched between 
two leads. Electrons are allowed to tunnel from the left lead (emitter) to the right lead (collector) 
under the influence of an optical pump. We start with the following Hamiltonian 

H = H d + H l + H dtl + H dte + H T , (1) 

where the first term describes electrons in the InAs/GaAs QD . We assume that the quantum con- 
finement effect is strong for the small QD considered here. Therefore, the energy spacings between 
the ground state and the first excited state for electrons and holes, AE e and AE h , are much larger 
than thermal energy, ksT, where fee and T denote the Boltzmann constant and temperature. Only 
the ground state levels for electrons and holes, E e and E h , are considered in H d . The second term 
describes the kinetic energies of free electrons in the electrodes, where the correlation effects among 
electrons is ignored. Note that in the current setup, the gate electrode does not provide any electrons, 
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but merely controls the energy levels of the QD. The third term describes the coupling between the 
QD and the leads. The fourth term describes the interband optical pumping with a frequency of uj , 
which is in resonance with the energy difference between an electron level in the wetting layer and the 
hole ground state level. Due to the large strain-induced splitting between the heavy-hole and light- 
hole band for typical QDs, we only have to consider the heavy hole band (with J z = ±3/2) and ignore 
its coupling with light-hole band caused by the QD potential. Because the effect of inter-particle 
Coulomb interactions is significant in small semiconductor QDs, we take into account the electron 
Coulomb interactions and electron-hole Coulomb interactions in the last term. 

Once the Hamiltonian is constructed, the tunnelling current of SET can be calculated via the 
Keldysh Green's function method. 6 We obtain the tunnelling current through a single dot 

J = ~TT 1 1^' ~ ^ ~ ~ ^» r^lt) /mG - (e) - (2) 

Eq. (2) is still valid for the SET under optical pumping provided that the condition >> R e h 

is satisfied, where R e h is the electron- hole recombination rate. /z,(e) and /^(e) are the Fermi distri- 
bution function for the source and drain electrodes, respectively. The chemical potential difference 
between these two electrodes is related to the applied bias via /il — = eV a . r^(e) and r#(e) 
denote the tunnelling rates from the QD to the left (source) and right (drain) electrodes, respectively. 
For simplicity, these tunnelling rates will be assumed energy and bias-independent. Therefore, the 
calculation of tunnelling current is entirely determined by the spectral function A = ImG r e a {e), which 
is the imaginary part of the retarded Green's function G r e a {e). 

The expression of retarded Green's function ,Gl a (e), can be obtained by the equation of motion 
of G r ea (t) = —i6(t){{d eja (t),dl a (0)}), where 9(t) is a step function, the curly brackets denote the 
anti-commutator, and the bracket (...) represents the thermal average. After some algebras 5 , the 
retarded Green's function of Eq. (2) 

GIM = (1 - Ne,-,){ 1 " K ' + " k w -} X nK ° nh -° P) 

+ n^g + n h _ a - 2n h>a n h _ a n h ^n h _ a 

e - E e + U eh + i^ e-E e + 2U eh + i^f 

, i» j A - ( n h,a + rih-a) + n h , a nh,-a 



nh,a + n h _ a - 2n ht(T n h - a ^ n hyCr n h _ a 

e - E e - U e + U eh + i^f e-E e -U e + 2U eh + i\ 



In Eq. (3), T e is the electron tunnelling rate T e = Tl + Tr. It is worth noting that the electron-hole 
recombination effect have not been directly included into Gg jCr (e), because we assumed that T e >> R eh . 
The typical value of R e h for InAs/GaAs QDs is ~ l/ns. The electron occupation number of the QD 
can be solved in a self-consistently via the relation 

M f de T L f L (e)+T R f R (e) T 

J TT L L + L R 

N e>a is limited to the region < N e>a < 1. Eq. (4) indicates that the electron occupation numbers of 
the QD, N e _ a and N et<7 , are primarily determined by the tunnelling process. To obtain the electron 
and hole occupation numbers (n e) _ CT = n e (7 and rih~ a = n^a) arised from the optical pumping, we 
solve the rate equations and obtain 

_ _ _ 7e,c-^e,fc /p-N. 

n e - n e ,_ CT - n e , CT - ^ ^ + R ^ + ^ , ^ 

and 

n h = n h - a = n h , a = —— — — — — , (6) 

lh, c N h , k + R eh {n e + N e ) + T h 

where r y e (h),c and N e (h) t k denote the captured rate for electrons (holes) from the wetting layer to the 
QD and the occupation number of electrons (holes) in the wetting layer. Here, we assume that N e ^ ;k 
is in proportion to the intensity of excitation power, p exc - T/i denotes the nonradiative recombination 
rate for holes in the QD. 

According to Eq. (3), particle Coulomb interactions will significantly affect the tunnelling current of 
SET. To illustrate this effect, we apply our theory to a self-assembled InAs/GaAs QD with pyramidal 
shape. First, we calculate the inter-particle Coulomb interactions using a simple but realistic effective- 
mass model. The electron (hole) in the QD is described by the equation 
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#e(h)^e(h)( r )> ( 7 ) 
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where m* e {p,z) (a scalar) denotes the position-dependent electron effective mass, which has m* eG = 
0.067m e for GaAs and m* eI = 0.04m e for InAs QD. m* h (p, z) denotes the position-dependent effective 
mass tensor for the hole. It is a fairly good approximation to describe m* h {p, z) in InAs/GaAs QD as 
a diagonal tensor with the x and y components given by m* t ~ l = (71 + ^2)jm e and the z component 
given by mp 1 = (7! — 27 2 )/m e . 7! and 72 are the Luttinger parameters. Vq D {p,z) (Vn D (p,z)) is 
approximated by a constant potential in the InAs region with value determined by the conduction- 
band (valence-band) offset and the deformation potential shift caused by the biaxial strain in the QD. 
These values have been determined by comparison with results obtained from a microscopic model 
calculation 7 and we have Vq D = — 0.5eV and Vq D = — 0.32eV. The eFz term in Eq.(7) arises from 
the applied voltage, where F denotes the strength of the electric field. Using the eigenfunctions of 
Eq. (7), we calculate the inter-particle Coulomb interactions via 

V iJ= j ir Ji* /^l^f\ (8) 

J J Co| r l — r 2| 

where = e, h. rij(ri) denotes the charge density. e is the static dielectric constant of InAs. The 
Coulomb energies are different in different exciton complexes, but the difference is small. 5 Therefore, 
only the direct Coulomb interactions have been taken into account in this study. 

For the purpose of constructing the approximate wave functions, we place the system in a large 
confining cubic box with length L. Here we adopt L = 40nm. The wave functions are expanded in a 
set of basis functions, which are chosen as sine waves 

i>nim(p,4>,z) = -^=sin(kix)sin(k m y)sin(k n z), (9) 
V L 6 

where k n = nn/L,k m = mir/L,ke = £tt/L. n, m and £ are positive integers. The expression of the 
matrix elements of the Hamiltonian of Eq. (7) can be readily obtained. In our calculation n = 20, 
m — 10 and £ = 10 are used to diagonalize the Hamiltonian of Eq. (7). Fig. 1 shows the inter- 
particle interactions as functions of QD size. The ratio of height and base length is h/b = 1/4, while 
h varies from 2.5nm to Q.hnm. The strengths of Coulomb interactions are inversely proportional 
to the QD size. However as the QD size decreases below a threshold value (around b = 12nra), 
U e is significantly reduced due to the leak out of electron density for small QDs. These Coulomb 
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interactions approach approximately the same value in the large QD limit. This indicates similar 
degree of localization for electron and hole in large QDs. We also note that U e h is smaller than U e 
in the large QD. This is due to the fact that in large QDs the degree of localization for the hole 
becomes similar to that for electron, while the anisotropic nature of hole wave function reduces U e h- 
The repulsive Coulomb interactions, U e and Uh, are the origin of Coulomb blockade for electrons and 
holes, respectively. The attractive Coulomb interaction U eh gives rise to the binding of the exciton. 
To study the behavior of tunnelling current, we consider a particular pyramidal InAs/GaAs QD 
with base length b = 13nm and height h = 3.5nm. The other relevant parameters for this QD are 
E efi = -O.UeV, Ehfl = -0.125eV,U e = W.lmeV, U eh = !6.7meV and U h = 18.5meV. 

Now, we perform detailed numerical calculation of the tunnelling current. For simplicity, we 
assume that the tunnelling rate T L = T R = 0.5meV is bias-independent. We apply a bias voltage 
V a across the source-drain and V g across the gate-drain. The QD electron and hole energy levels, 
E e and E h) will be changed to E e + aeV a — f3eV g and E h + aeV a — f3eV g) where a and (3 are the 
modulation factors. In our calculation, we assume a = 0.5 and j3 = 0.7, which can also be determined 
by experiments 9 . Meanwhile the chemical potentials of the electrodes with Fermi energy Ep = 60meV 
(relative to the conduction band minimum in the leads) are assumed to be 70meV below the energy 
level of E e at zero bias. Parameters T h = 0.2meV and R e h = lOfieV are adopted. 

Applying Eqs. (2), (4) and (6), we solve for the electron occupation number N e = N ej(T = N e _ a 
and tunnelling current J. Fig. 2 shows the calculated results for N e and J as functions of gate 
voltage with and without the photon-excitation power at zero temperature and V a = 2mV. Solid 
line and dashed line correspond to / = (no pump) and / = 0.9 (with pump), respectively. We 
have defined a dimensionless quantity, / = ^h^Nh^/^h, which is proportional to the pump power. 
The electron occupation number displays several plateaus, while the tunnelling current displays an 
oscillatory behavior. We label four critical voltages (from V g \ to V gi ) to indicate the resonance energies 
of retarded Green's function. We see that the photon-excitation leads to additional two current peaks 
below the voltage V g3 , which is caused by the electron tunnelling assisted by the presence of a hole in 
the QD. This interesting phenomenon was observed by Fujiwara et al. in an SET composed of one 
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silicon (Si) QD and three electrodes 9 . The behavior of the photo- induced tunnelling current can be 
understood by the analysis of the poles of retarded Green's function of Eq. (3); the first peak of dashed 
line corresponds to the tunnelling current through the energy level at e = E e — 2U e h (corresponding to 
a positive trion state) . The second peak is caused by a pair of poles at e = E e — U eh (the exciton state) 
and e = E e + U e — 2U eh (the biexciton state). Since the magnitude of U e is very close to that of U e h, 
these two poles almost merge together. The third peak is caused by another pair of poles at e = E e 
(the single-electron state) and e — E e + U e — U eh (the negative trion state). The last peak locating 
near V g = 123mV is due to the tunnelling current through the energy level at the pole e = E e + U e 
(the two-electron state). The gate voltage difference AV^i = V g 2 — V g i (^43 = V g 4 — V g z) determine 
the strength of electron-hole interaction /5AV^ 2 i = U eh ((3V g43 = U e ). Once U eh is determined, we 
obtain the binding energy of exciton complexes. 

In this study we have used the tunnelling current of an SET under optical pumping to determine 
the electron-hole interaction U e h, which can be regarded as the binding energy of exciton. Although 
we used InAs/GaAs SET as an example, this idea can also be applied to Si/Si0 2 SET system 9 . 

ACKNOWLEDGMENTS 

This work was supported by National Science Council of Republic of China under Contract Nos. 
NSC 93-2215-E-008-014 and NSC-93-2120-M-008-002 



7 



REFERENCES 

L. Landin, M. S. Miller, M. E. Pistol, C. E. Pryor, and L. Samuelson, Science, 280, 262 (1998). 

P. Michler, A. Kiraz, C. Becher, W. V. Schoenfeld, P. M. Petroff, L. D. Zhang, E. Hu and A. 
Imamoglu: Science 290, 2282 (2000). 

Z. Yuan, B. E. Kardynal, R. M. Stevenson, A. J. Shield, C. J. Lobo, K. Cooper, N. S. Beattlie, D. 
A. Ritchie and M. Pepper: Science 295, 102 (2002). 

C. Santori, M. Pelton, G. Solomon, Y. Dale, and Y. Yamamoto, Phys. Rev. Lett. 86, 1502 (2001). 
David. M. T. Kuo and Y. C. Chang, Phys. Rev. B 69, 041306 (2004). 

H. Haug and A. P. Jauho, Quantum Kinetics in Transport and Optics o Semiconductor (Springer, 
Heidelberg, 1996). 

O. Stier, M. Grundmann, and D. Bimberg, Phys. Rev. B 59, 5688 (1999). 
J. Weis, R. Haug, K. Klitzing, and K. Ploog, Phys. Rev. B 46, 12837 (1992). 
A. Fujiwara, Y. Takahashi and K. Murase, Phys. Rev. Lett. 78 , 1532 (1997). 



8 



Figure Captions 

Fig. 1: Intralevel Coulomb interactions U e and U h and interlevel Coulomb interaction U e h as a 
function of the QD base length b. 

Fig. 2: Electron occupation number N e and tunnelling current as functions of gate voltage at zero 
temperature for various strengths of optical excitation. Current density is in units of Jo = 2e x meV/ h. 
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